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We have considered a mechanism for inducing a time-reversal violating electric dipole moment 
(EDM) in atoms through the interaction of a nuclear EDM djv with the hyperfine interaction, the 
"magnetic moment effect". We have derived the operator for this interaction and presented ana- 
lytical formulas for the matrix elements between atomic states. Induced EDMs in the diamagnetic 
atoms ^^^Xe, '^^'^Yb, ^^^Hg, ^^^Rn, and '^^^Ra have been calculated numerically. From the experi- 
mental limits on the atomic EDMs of ^^^Xe and ^^^Hg, we have placed the following constraints on 
the nuclear EDMs, |djv(^^^Xe)| < 1.1 x 10"^^ |e| cm and |djv(iS''Hg)| < 2.8 x IQ-'^* |e|cm. 
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I. INTRODUCTION 

Since the discovery of CP-violation in the decay of 
long-lived mesons, there have been many efforts to 
observe CP or time-reversal (T) violation in other sys- 
tems. The latter is equivalent to CP-violation assuming 
the validity of the CPT theorem. In particular, the pos- 
sible existence of a permanent electric dipole moment 
(EDM) of a particle would imply the violation of both 
parity and time- reversal invariance; see, e.g., [1]. The 
EDMs of the particles predicted by the standard model 
are too small to be detected at the present level of exper- 
imental accuracy. However, different extensions of the 
standard model (such as supersymmetry) predict much 
larger EDMs of the particles that could, in principle, be 
found using modern experimental techniques, see, e.g.. 

Recently, Griffith et al. reported a 7-fold im- 
provement on the limit of the atomic EDM of ^^^Hg, 
\d{^^^Rg)\ < 3.1 X 10-29 |e| cm. This is the most strin- 
gent limit on an atomic EDM. In order to interpret the 
measurement in terms of fundamental P,T- violating pa- 
rameters, atomic calculations are required. In our recent 
paper Q (see also ^-i^) we calculated atomic electric 
dipole moments induced by the nuclear Schiff moment, 
the P,T-odd electron-nucleon interaction, and the elec- 
tron electric dipole moment and placed limits on the cou- 
pling constants of the P,T-odd interactions from the new 
Hg result. 

In this work we consider one more P,T-odd interaction 
which gives rise to an atomic EDM in the second order of 
perturbation theory. This interaction was first discussed 
by Schiff in Ref. Q. He showed that if a nucleus has a 
permanent EDM dpf, an atomic EDM may be induced 
due to the interaction of djv with the magnetic field cre- 
ated by the electrons at the nucleus. 

It is worth noting that the existence of a non-zero djv 
alone is insufficient for producing observable EDM effects 



in neutral atoms due to electronic screening of an applied 
electric field at the nucleus @ . This screening is exact for 
the case of a point-like nucleus experiencing electrostatic 
forces. One can circumvent this screening by accounting 
for the finite size of the nucleus or the hyperfine interac- 
tion. The former leads to the appearance of the nuclear 
Schiff moment; the latter "magnetic moment effect" is 
the subject of this work. 

In his work 9], Schiff evaluated the EDMs of H and 
He induced by the magnetic moment effect. Hinds and 
Sandars later calculated the effect in TIF (T^]. In the 
current work we present a general analysis of the effect for 
atoms and perform calculations for diamagnetic atoms of 
experimental interest. 

The volume effect (nuclear Schiff moment) is gener- 
ally considered to be the dominant mechanism induc- 
ing EDMs in heavy diamagnetic atoms with nuclear spin 
/ = 1/2. However, it has been discovered recently that 
the nuclear Schiff moment is very sensitive to many- 
body corrections, see, e.g., [11,] and references therein. 
These corrections suppress the bare values for the Schiff 
moments for all considered nuclei. Other than a gen- 
eral suppression, there is yet no agreement between the 
many-body approaches. The nuclear EDM contribution 
to the atomic EDM (through the magnetic moment ef- 
fect) should therefore not be disregarded before the nu- 
clear many-body problem is well-understood and specific 
CP-violating models considered. 

Further motivation for this work comes from the grow- 
ing interest in measuring nuclear EDMs in ion storage 
rings [T2| - [l5j . In ions the nuclear EDM is not screened 
and can be measured directly. The relations obtained 
in this work between the nuclear and atomic EDMs en- 
able one to place limits on nuclear EDMs from neutral 
atom measurements. These limits may be considered as 
an accuracy benchmark for proposed nuclear EDM mea- 
surements with ions. 

The paper is organized as follows. In sections |TT] and 
mil we derive analytical expressions for the P,T-odd op- 
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erator and the matrix elements of this operator between 
atomic states. In Sec. II VI we present the equation for the 
atomic EDM and discuss different contributions to 
d^. In Sec.|V]we obtain simple analytical formulas that 
can be used for an estimate of the atomic EDM (and this 
is compared to the atomic EDM induced by the nuclear 
Schiff moment), we describe our numerical method for 
EDM calculations of diamagnetic atoms, and we present 
results for ^^^Xe, i^^Yb, ^^^Hg, ^URn, and ^^'^Ra. Sec- 
tion |Vl] contains concluding remarks. 



II. GENERAL FORMALISM 

Let us assume that the nucleus has a P,T-odd EDM, 
dN = (dw) = dN^/I- Here (dAr) denotes the expectation 
value of the nuclear electric dipole moment with the exact 
nuclear ground state wave function Ig*] . 

The P,T-odd operator corresponding to the magnetic 
moment effect can be written as [9| 



U = -t[Q,HM], 



(1) 



where square brackets in Eq. ([T]) denote a commutator. 
The operator Q is determined as (if not stated otherwise 
we use atomic units h — m = \e\ = I) 



Q 



(2) 



fc=l 



where Z is the nuclear charge and pn and Pe ' are mo- 
mentum operators for the nucleus and electrons, respec- 
tively. Hm is the operator for the hyperfine interaction 
(HFI) which may be interpreted as the interaction of the 
nuclear magnetic moment fij^ with the magnetic field of 
the electrons. It can be represented by the sum of the 
single-electron operators 



Vi X OLi 



(3) 



are the Dirac matrices. fJ-N — I^n^/I 



where ri> = max(ri,i?) with R the nuclear radius, and 

(T 

a 

and fJ-N ~ /^/(2mpc), where fi is the nuclear magnetic 
moment, nip is the nucleon mass and the speed of light 
c= l/a« 137. 

A note on the origin of the P, T-odd operator Eq. ([Ij . 
Schiff showed (we limit the explanation to the case of an 
atom with nuclear EDM) that when an atom is acted 
upon by electrostatic forces only and the nuclear EDM 
and charge distributions are the same, the full Hamilto- 
nian may be expressed as H — Hq + i[Q, Hq], where the 
full Hamiltonian H includes the nuclear EDM and Hq 
does not. The eigenvalues of H are therefore the same as 
those of Hq to first order (the expectation value of the 
commutator containing Hq is zero) and do not contain 
the nuclear EDM, that is, there is no linear Stark shift 



and no observable EDM of the atom. When the hyper- 
fine interaction is taken into account, the full Hamilto- 
nian may be written as H' = Hq + i[Q, Hq] — i[Q, Hm], 
where the prime denotes that the magnetic interaction 
is included. Only the commutator containing Hm may 
lead to observable EDM effects. We refer the reader to 
Schiff's landmark work Ref. [sl for details. 

In Eq. ([2]) we assumed that the center of mass is at 
rest. As a result, the momentum operator of the nucleus 
Pn can be replaced by the sum of the electronic momenta 

(k) 

Pe with the opposite sign. In the following we will deal 
with the electronic momenta only and omit the subscript 
e, i.e., we denote p = Pe. Then we obtain for the operator 
U: 



Z r 



ZP 



k=l 



rfc X Oik 



'k> 



(4) 



fc=i 



We will consider the single-electron operator Uk (omit- 
ting for brevity the index k) since, as seen from Eq. (j4|. 
a generalization to the case of a many-electron atom is 
straightforward. In addition, we restrict ourselves to con- 
sideration of the most interesting case of nuclear spin 
/ = 1/2. Then for the Cartesian components m and i of 
the nuclear spin I we have 



(5) 



Substituting Eq. ([5]) in Eq. (|4]) and taking into account 
that the terms ^ Smi will be canceled out, after simple 
transformations we obtain 



U = 







f CTArr 1 


acTN |p, ^1 








f 





(6) 



Here (Tn — 21 and {...}^_ is an anticommutator. Now we 
take into account that 

d n X L 

P = -»V = -zn- , 

or r 

with L = r X p the orbital momentum operator. 
Hence 



'> J + 
and finally we obtain 



d r 



U = - 



Z 



d r 
dr ' r5 



. ap 



(7) 



If the typical distances of interest are small one can 
neglect the eigenvalue and the mass of the electron in 
comparison with the electrostatic potential. Hinds and 
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Sandars showed in Ref. [l^] that in this approximation 
the operator U can be written in a more simple form: 



2 dNl^^N a X L 
^ = 7^ — ""w 5 — ■ 



(8) 



In the following sections we will discuss the difference 
between the two forms of the operator U given by Eqs. 
dZD and ©. 



IV. ATOMIC EDM 

In this section we present our calculations of EDMs 
induced by the operator U for several diamagnetic atoms 
of experimental interest having I = 1/2. The EDM = 
(TN of an atom in the state |0) in the second order of 
perturbation theory is given by 

(0|D|if) {K\U\Q) 



K 



Eo — Ek 



(12) 



III. ELECTRONIC MATRIX ELEMENTS 

We use the following form for the electronic wave func- 
tions 



\nKm) = 



ig{r)n 



where I = 2j — I. Using the expression for the operator 
U given by Eq. ([7]) we can derive the electronic matrix 
elements (MEs). Finally we find (see the Appendix for 
details of the derivation) 



{n' x'm! \ U\ n^cm) = — 



Z 



■ (Tpf{>c m \n\ xm) 



{fg{x^x' + l) + fg'{x'^x+l)} 



+ (e'-e) {x-x'){ff-g'g)]—dr. 



(9) 



If we factorize (e' — e) in the second term of Eq. ^ , 
we can rewrite it as 



(e — £ ) — - — (Tn (n X m. 



Then Eq. (|9]) can be represented by 

dNfJ-N 



nxmj 



(10) 



{n' x'm! \U\ nxm) 



criq[{x' m' \Ta.\ xm) x 



{x+x') / {f'g{x- x' + 1) + fg' {x' - x+l)} 



+ (e' — e) ( n' x'm 



r X cr 




y,3 


nxm^ 



< > 

(11) 



As seen from Eqs. ^ and (|lip . the first term in these 
equations disappears if x' + x = 0. It happens for the 
states with I' = I ±1 and j' ~ j. As a result, the MEs 
{—xm' \U\ xm) turn out to be proportional to (e' — e). 
For heavy atoms these MEs are small. In particular, it 
means that the MEs (s|J7|pi/2) contribute less to the 
atomic EDM than the MEs (s If/] P3/2). If we neglect 
the term (e' - e) in Eqs. © and (HH), then Eqs. © 
and ([5]) lead to the same formula for the MEs. 



where D = — r is the electric dipole operator and the 
summation goes over all intermediate states \K) allowed 
by the selection rules. 

As follows from the consideration given in the preced- 
ing section, it is convenient to present the operator U as 
a product of two operators, one relating to the electronic 
part and the other relating to the nuclear part, i.e. 



U — Uci <TAr, 



(13) 



where Ud describes the electronic part of the operator 
U. 

Applying the Wigner-Eckart theorem, and summing 
over magnetic quantum numbers to of the initial and final 
states in Eq. (fT^ . we obtain for d. 



N . 



^ _2^{Q\\r\\K) mU.xWK) 

"at ~ o / . 



K 



Ek — Eq 



(14) 



To carry out calculations of atomic EDMs for atoms 
with closed shells, it is convenient to rewrite Eqs. (IT^ and 
([T^ in terms of single-electron wave functions. Then the 
corresponding expressions for and are given by 



2E 



at '-^^^'-^ '-''at 

(c|r|n) {n\U\c) 



and 



(15) 



(16) 



r^dr where the indices c and n relate to the single-electron 
core and virtual orbitals and £c(n) £^re the single-electron 
core (virtual) energies, correspondingly. 

Note that if we substitute the second term of Eq. (fTT|) 
into Eq. then after cancelation of {£n~£c) in 

the numerator and denominator and applying closure 
~ 1' obtain a term which is proportional 

to 



E 



r [(r X cr) <T 



It can be readily shown that if we apply again the 
Wigner-Eckart theorem and sum up over the projections 
TOc of the total angular momenta jc, this term goes to 
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zero and therefore does not contribute to d^. In partic- 
ular, it means that the second term in Eq. (jlip does not 
contribute to the atomic EDM for closed-sheU atoms. 

Keeping in mind that we intend to carry out calcu- 
lations of EDMs for atoms with closed shells, we can 
neglect the term ^ (e' — e) in Eq. (fTT]) . Accounting for 
Eq. (jl3|) . we obtain the following expression for the re- 
duced ME of the operator Ud'. 



(nV'l |(7ei| \nK) 



— I 



X / {S'g{>c->c' + 1) 
Jo 

flT 



where the reduced ME {x'\\n\\>c) is given by 



(17) 



(V| |n|| x) = {-ly'+'f^ V(2/ + l)(2j + l) 

>^(-{/2 1/2 + ^ + (18) 



with 



1, if a; is even 
0, if X is odd 



In Eq. ([TTI) we also took into account that the integral 
inside the nucleus is very small and replaced r> by r. As 
we mentioned above, the MEs of the operator L^ci turn 
to zero \i K + = Q. It means that MEs between states 
with the same total angular momentum though different 
parity like (s||J7ci|bi/2), (P3/2 1 1 1^0111^3/2), etc. are equal 
to zero. 



V. METHOD OF CALCULATION AND 
RESULTS 

A. Analytical estimates 

In this section we derive the analytical expression for 
the matrix element Eq. p7|) . Outside the nucleus the 
wave functions and gn>c can be represented by [l7j 



fn 



1 



(7 + x)J2-i{x) - -J2',-l{x) 
ZaJ2-i{x), 



(19) 



where x = \/8Zr, 7 = \/ — Z'^c? ^ and Ju{x) are Bessel 
functions. 

Using these expressions for the radial wave functions 
we obtain for the radial integral 



fn' yc' Qr\ 



dr 



Z^a 



3/2 



\R 



■M- 



(20) 



Here 



A = 



-96 



1 



1 



A, 



(21) 



the relativistic enhancement factor Rm is given by [I| 
1 r(7' + 7-2) 



Rm — 



Aj,j r(7' - 7 + 3)r(7 - 7' + 3)r(7' + 7 + 3) ' 

is deter- 



T{j3) are the 7-functions, and the factor Ajij 
mined as follows 



(/+J-2)! 



(/-J + 2)!(j-j' + 2)!(f +J+3)!- 



(22) 



Let us consider an important particular case of the ME 



|C/cl||P; 



'3/2/ 



In this case x' — —1 and x — —2. Then 



the first term under the integral in Eq. (|17p turns to zero 
while for the second term in the integral we find, using 
Eqs. (|20 l) and ([2T 1) . 



dr 



Z^a 



,3/2 



R 



M- 



For the matrix element we therefore have 



l's||C/el||"-P3/2) ~ dNfJ-N- 



2y/3Za 



Rm- 



(23) 



(24) 



This analytical expression demonstrates how the atomic 
EDM depends on Z and nuclear parameters. The en- 
hancement factors Rm for the medium atoms are close 
to unity. They grow with increasing Z approaching the 
value 2.2 for Ra. The values of the enhancement factors 
for the diamagnetic atoms considered in this work are 
listed in Table HI 

A simple analytical estimate gives the following result 
for the atomic EDM 



±10-' fiZRMd 



(25) 



where the upper sign {+) relates to the divalent atoms 
and the lower sign (-) relates to the noble gases. This 
formula gives values in reasonable agreement (within a 
factor of ~ 2) with the many-body Dirac-Hartree-Fock 
(DHF) results presented in the following section. 

It is instructive to compare this estimate with a similar 
one for the atomic EDM induced by the nuclear Schiff 
moment S, 

4 « T 10-22^2 (i?i/2 + 2i?3/2) ^/(|e| fm^) |e| cm , (26) 

where R1/2 and R3/2 are relativistic enhancement fac- 
tors corresponding to s — Pi/2 and s — P3/2 weak matrix 
elements, respectively; see Ref. [l^ for these factors. 
Again, the upper sign is for the divalent atoms and the 
lower sign is for noble gases. The values obtained from 
this simple formula also reasonably agree (within a fac- 
tor of 2) with the many-body DHF results obtained in 
Refs. 
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TABLE I. The nuclear charges Z and the relativistic enhance- TABLE IL The values of in units (10~^ djv) obtained in 
ment factors Rm- the DHF and RPA approximations. 





i29Xe 






2iiRn 


Ra 




^^^Xe54 






'"Rnse 


'"''Rags 


z 


54 


70 


80 


86 


88 


DHF 


4.4 


2.6 


3.2 


-9.1 


-9.5 


Rm 


1.28 


1.56 


1.83 


2.06 


2.15 


RPA 


5.8 


11 


11 


-13 


-33 



The atomic EDM induced by the Schiff moment ben- 
efits from an extra Z dependence which becomes very 
important in heavy atoms 17[. In Sec IV CI we consider 



a specific mechanism for inducing the Schiff and nuclear 
dipole moments and compare the sizes of the induced 
atomic EDMs. 



B. Numerical method of calculation and results 

Here we describe the simple numerical methods we use 
for calculations of atomic EDMs for closed-shell atoms. 
At the first stage we solve DHF equations in the 
approximation. This means that we include all electrons 
forming the ground state of the atom in a self-consistency 
procedure 



(27) 



Here iJo is the relativistic Hartree-Fock Hamiltonian and 
"^c are single-electron wave functions of the core. 

To take into account polarization of the atomic core 
by external fields (the electric dipole field or the P,T-odd 
field), we solve the random phase approximation (RPA) 
equations 



(i/o 



=)(5?/'c 



(28) 



where F is the operator of the external field and SV^ is 
the correction to the self-consistent potential due to the 
effect of the external field. The RPA equations are 
solved self-consistently for all states in the core. 

In implementing the DHF and RPA procedures for cal- 
culations of the atomic EDMs, we have used two equiv- 
alent approaches. The first involves the construction of 
virtual orbitals and summation over states. The second 
involves the direct solution of the perturbed orbital Jt/'c 
on the grid. 

In the former method the virtual orbitals are con- 
structed by multiplication of the previous orbital of the 
same partial wave to a smooth function of r with subse- 
quent orthogonalization of this orbital to the rest of the 
orbitals. This method was described in detail in Refs. fl^ 
and |2a]. 

In the latter method, implemented in Ref. for cal- 
culation of the Schiff moments, instead of direct summa- 
tion over virtual states in Eq. (1151) . we evaluate = 

2Ec(c| - where |^c^) = "Sri") ^ ^o^^" 

tion of Eq. (^5)) for the P,T-odd field. In taking into ac- 
count core polarization by the fields, the correction goes 



to one field and not the other to avoid double counting. 
See Ref. [7| for details. 

The results of numerical calculations carried out using 
the DHF and RPA methods for the considered diamag- 
netic atoms are presented in Table [Hi It is seen that 
inclusion of the RPA corrections increases the size of the 
atomic EDM. For the noble gases (Xe and Rn) the RPA 
corrections contribute at the level of 30-40%, while for 
atomic Hg, Yb, and Ra, which have two s electrons above 
closed shells, the RPA corrections are much larger. In 
fact, they increase the EDMs of these atoms several times 
compared to the DHF values. The reason for this increase 
is that the two s electrons are loosely bound and can 
be easily excited. As a result, account of higher orders 
of perturbation theory (like the RPA corrections) leads 
to a significant change in the "bare" results obtained in 
the DHF approximation. This result is similar to that 
obtained for other contributions to the P,T-odd atomic 
EDM discussed in Refs. 0, Q for the same atoms. 

We note that in our previous works [1, 0|, a check 
on the RPA results for Yb, Hg, and Ra was carried 
out by performing more sophisticated calculations in the 
V^^"^ approximation. In this approach, correlations 
between the two valence electrons and the core were 
taken into account using many-body perturbation theory 
(MBPT) while correlations between the valence electrons 
were accounted for using the configuration interaction 
(CI) method; see Refs. lluil regarding the CI-I-MBPT 
method. It was found in [3, Q] that the two approaches 
{y^ and y^-2) yield results for the EDMs that differ 
by less than 20% for the various P,T-odd mechanisms. 
The good agreement of the results obtained using two 
very different approaches is a strong argument in favour 
of the stability of the RPA results. Because the operator 
considered in this work is of similar form to the other 
P, T-odd operators, we expect that our RPA results for 
the atomic EDMs are accurate to about 20%. 

Using the experimental limits on the P,T-odd atomic 
electric dipole moments of ^^^Xe [23j 



d(129Xe) = (0.7±3.3stat±0.1syst 

M(^^^Xe)| < 6.6 X 10"^^ 



X lO"^'^ e cm 



e cm 



(29) 



and is^Hg [ 



(0.49 ± 1.29,tat ± 0.76,y,t) x IQ-^^ e cm 
[g)| < 3.1 X 10"^^ |e|cm (30) 



we are able to place constraints on the nuclear EDMs 
dAr(i29Xe) and dt^Q^^^^g). Using Eqs. ([Ml), O, and 
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the results presented in Table [TT] we obtain 
|dw(i29xe)| < 1.1 X 10-21 |e|cm, 
\dNC^^ilg)\ < 2.8 X 10" 



e cm. 



(31) 



We have not included the atomic theory error in these 
limits. 



C. Contributions to the nuclear EDM 

For spherical nuclei with spin determined by a single 
unpaired nucleon, there are several main terms that con- 
tribute to the nuclear EDM. One of them is character- 
ized by the P,T-odd nucleon-nucleon interaction while 
the other ones are the contributions from the EDMs of 
the neutron (c?„) and the proton (dp). Thus, we can write 
dN as 



dN = 



tjd„ 



Qpdp , 



(32) 



where we denote by d]^ the contribution from the P,T- 
odd nucleon-nucleon interaction; tj = 1 ior I — Ij + 1/2 
and ti = -I /{I + 1) for I = Ij - 1/2 (with // being 
the orbital momentum of the unpaired nucleon) and the 
coefficient a„ is numerically close to 0.1. 

In Refs. [ll, m it was shown that the nuclear EDM 
induced by the P,T-odd nucleon-nucleon interaction can 
exceed the nucleon EDM by more than two orders of 
magnitude. If we neglect the terms proportional to dn 
and dp in Eq. (f32|) and express through the T-odd 
nucleon-nucleon coupling constant rj, we obtain (23 | 



I TV 



4 X 10"^^ 



(33) 



where q — and 1 for an outer neutron and proton, 
respectively. Taking into account that for all atoms con- 
sidered q = 0, we obtain for d^ (in units of |e| cm) 



cJat « — 2 X 10 |e| cm. 



Z 



(34) 



Interactions of the outer neutron with both protons and 
neutrons of the core contribute to the nuclear EDM; the 
coupling constant rj = ^rynp + ^f?nn. Using Eq. ([M)) 
and the results for given in Table [Tl] we can express 
the values of the atomic EDMs through rj. These results 
for i^^Xe, ^^^Yb, ^^^Hg, and ^^^Rn are presented in Ta- 
ble |TTI1 

To get an idea of the relative size of the induced atomic 
EDMs compared to those induced by the Schiff mo- 
ment, we consider the Schiff moment produced by the 
nucleon-nucleon interaction. For all considered atoms, 
the nuclei have an outer neutron. This means that there 
is no direct contribution to the Schiff moment. The 
protons need to be excited to distort the charge den- 
sity and create a P,T-odd charge distribution. A sim- 
ple numerical calculation in the Woods-Saxon poten- 
tial with spin-orbit interaction included gives for ^^^Hg, 



TABLE III. is the spin and parity of the nuclear ground 
state and /i is the magnetic moment expressed in nuclear mag- 
netons 29]. The values of in units (10"^^ |e| cm) are 
obtained in the RPA approximation. 





F 




ti 




i^«Xe 


1/2+ 


-0.7780 


1 


-5.1 


"^Yb 


1/2- 


0.4919 


-1/3 


3.1 


199Hg 


1/2- 


0.5059 


-1/3 


3.1 


2iiRn 


1/2- 


0.60 


-1/3 


-3.7 


225Ra 


1/2+ 


-0.734 


1 





S'(iayHg) = -1.4x 10-Vp|e| fm-' [2|. Using this resuh, 
the calculation for the atomic EDM induced by the Schiff 
moment df^ = -2.8 x 10-^^ S/{\e\ fm^) |e| cm [IQ, and 
the result for the atomic EDM induced hy dN through 77 
presented in Table [1111 we obtain 



dgr^ng) 



df,(i99Hg) 



0.01 



(0.477np + 0.6r? 

n 

Vnp 



(35) 



(The ratio is larger for lighter atoms.) While this in- 
dicates that the contribution to the atomic EDM from 
the nuclear EDM is significantly smaller than that from 
the nuclear Schiff moment, we remind the reader that we 
have used a very simple model for the nucleus. Indeed, it 
is only recently that many-body calculations have been 
performed for the nuclear Schiff moment and these have 
demonstrated that many-body corrections are large and 
lead to a suppression of the bare results for all nuclei con- 
sidered, see, e.g., the most recent calculation [llj and ref- 
erences therein. The results of the different many-body 
approaches, however, are not in agreement. For exam- 
ple, in their RPA approach, Dmitriev and Sen'kov [2^ 
find a suppression of two orders of magnitude in the 
isoscalar channel of the P, T-odd pion-nucleon-nucleon in- 
teraction and a suppression of one order of magnitude 
in the isotensor channel for the case of ^^^Hg. In the 
fully self-consistent approach of Ban et al. [ll| applied 
to ^^^Hg, one order of magnitude suppression is seen in 
the isoscalar and isotensor channels, while instabilities 
are seen in the isovector channel, with results even vary- 
ing in sign. Therefore, until the many-body problem is 
well-understood, and specific CP-violation models con- 
sidered, the contribution to the atomic EDM from the 
nuclear dipole moment should not be dismissed. 

It is worth noting that there are nuclei with octupole 
deformation, such as, e.g., ^^^^Rn and ^^^^^^^^Ra. For these 
nuclei, the nuclear EDM cannot be approximated by the 
simple formula Eq. As shown in [26, 27], the P,T- 

odd nuclear forces lead to an enhanced collective dipole 
moment that can significantly exceed single-particle mo- 
ments. The enhancement of the nuclear EDM is also 
possible in nuclei with quadrupole deformations due to 
mixing of close opposite-parity levels [l^ HI] . 
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VI. CONCLUSION 

We have derived an expression for the P,T-odd op- 
erator produced by the interaction of a P,T-odd elec- 
tric dipole moment of the nucleus with the operator 
of the hyperfine interaction. We have presented sim- 
ple analytical formulas that can be used for an esti- 
mate of the EDM for different atoms. Using numeri- 
cal methods, we have found the contributions to atomic 
EDMs caused by this P,T-odd interaction for a num- 
ber of diamagnetic atoms. Using the experimental lim- 
its on the atomic electric dipole moments of ^^^Xe and 
i99jjg^ we constrain the EDMs of the nuclei of ^^^Xe 
and i99Hg to be \dN{^^^Xe)\ < 1.1 x 10-2i|e|cm and 
|d7v(^^^Hg)| < 2.8 X 10"^'* |e| cm, correspondingly. These 
limits give an accuracy benchmark for proposed measure- 
ments of nuclear EDMs in ion storage rings 

This work was supported by the Australian Research 
Council. We would like to thank G. Jones for help at the 
initial stage of the work. 



Appendix A 

For the matrix elements of the two operators enter- 
ing Eq. ([T]) (see the main text) we obtain after certain 
transformations 



Adding (|Aip and (|A2p , we find for the ME of the operator 
U: 



{n' x'm! \ U\ riKm) = — 



Z 



-CT/v [x m \n\ xm) 



(V - x) (g^ - /'^ - g'^ + f^] 
\ dr dr dr dr J 



x" 



ifg + fg') 



-dr. 



(A3) 



Using the Dirac equations for the radial wave functions 



df l + x', f , Za\ , 

^ = / +[e' + 171+ ]g\ 

dr r \ r J 

dg I — X f Za\ 

— = g - e - m H /, 

dr r \ r J 

df l + x ( Za\ 

— = / + £ + TO H U, 

dr r \ r J 

dg' 1 " x', ( . Za\ , , ^, 

-f = g' -ie' -m+ /'. (A4) 

dr r V T J 



d r 

r ri. 



(Tn (x'm' Inl xm) 



[x - x-l)[ g— - f -r 
dr dr 



+ -l)[g-^f — 



and 



crjvr 

— 3-,"P 



-^dr 

rt. 



(Al) 



cttv {x'm' |n| xm) 



^ifg + fg') 



dr dr dr dr 



-dr. 



(A2) 



we can write 



dr dr dr dr 



ifg'~fg)~{e'-s){f'f-g'g). (A5) 



Now substituting (|A5p to (|A3[) we finally arrive at Eq. © 
given in the main text. 
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